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Abstract 

Optimizing resonator shapes for maximizing the 
ratio of maximum to minimum gas pressure at an 
end of the resonator is investigated numerically. 
It is well known that the resonant frequencies 
and the nonlinear standing waveform in an 
acoustical resonator strongly depend on the 
resonator geometry. A quasi-Newton type 
scheme was used to find optimized axisymmetric 
resonator shapes achieving the maximum 
pressure compression ratio with an acceleration 
of constant amplitude. The acoustical field was 
solved using a one-dimensional model, and the 
resonance frequency shift and hysteresis effects 
were obtained through an automation scheme 
based on continuation method. Results are 
presented for optimizing three types of 
geometry: a cone, a horn-cone and a half cosine- 
shape. For each type, different optimized shapes 


were found when starting with different initial 
guesses. Further, the one-dimensional model was 
modified to study the effect of an axisymmetric 
central blockage on the nonlinear standing wave. 

INTRODUCTION 

The waveform of the standing wave in an oscillating 
closed cavity is strongly influenced by the geometry of 
the cavity (resonator). Lawrenson et ai 1 at MacroSonix 
Corp first exploited the dependence for developing 
high-amplitude acoustic pressures, referred to as 
resonant macrosonic synthesis (RMS). Peak acoustic 
pressures that exceed the ambient pressure three to four 
times and ratios of maximum to minimum pressures of 
27 were observed in RMS cavities. The size of the 
demonstrated overpressure reached the level that 
required by commercial applications such as acoustic 
pump or compressor. They considered these types of 
axisymmetric shape: cylinder, cone, hone-cone and 
bulb, and appeared to find that the hone-cone shape 
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generates higher overpressure given an input power. A 
companion paper by Ilinskii et al 2 developed a one- 
dimensional time domain model for studying the RMS 
numerically. Their results confirmed the nonlinear 
standing waveform and the related characteristics such 
as shape-induced resonance hardening and softening 
observed in the experiments by Lawrenson et al x Chun 
and Kim 2 numerically investigated the cosine shapes in 
addition to cylindrical and conical shapes using high- 
order finite -difference approximations. They claimed 
that the half cosine-shape is more suitable to induce 
high compression ratio than other shapes. To serve 
commercial needs, such as in acoustic gas compressor 
and acoustic liquid pump, the common goal in these 
studies is to find an optimized shape for generating 
higher overpressure. However, the optimization 
procedure and result have not yet been discussed. 

In this article, we introduce the numerical schemes 
and present the results in searching for the optimal 
shape parameters that give us the highest maximum-to- 
minimum pressure ratio in each of the following types 
of shapes: cone, horn-cone and half cosine-shape. In the 
second part, we study the effect of a cylindrical central 
blockage on the RMS. 


I. Governing Equations 


In this section, we describe a one-dimensional model 
for the acoustic wave field in the resonator that includes 
a central blockage. We develop our model following 
the work by Ilinksii et al. 2 for an axisymmetric 
resonator in absence of the central blockage. 

Consider the acoustic field in an oscillating 
resonator driven by an external force. As sketched in 
Fig. 1, the resonator is axisymmetric with the inner 
radius given by r=r(x ), 0<x<l , where x is the 
coordinate along the axis of symmetry. The central 
blockage is also axisymmetric with the same axis of 

symmetry x, whose shape is written as r 0 = /J,(x) . 


The density of the gas p, the velocity u and the pressure 
p satisfy the conservation of mass 


d P 


+• 


i 


- r Q 9X 


((r 2 -r 2 )pu) = 0. 


dt r 

and the conservation of momentum 


(1) 


du du _ 1 dp 

dt dx p dx 

(£ + 4?7/3) 8 f 


a(t) 


+ 


dx 


-r 0 2 dx 


((r 2 ~r 2 )u) 


A 2) 


where aft) is the acceleration of the resonator enforced 
by the external force; £, and rj are coefficients of 
viscosity. The state equation is specified by that of an 
ideal gas 


P = Po(P / PoY’ 0 ) 

where p Q and p 0 are the ambient pressure and density 

respectively. The no-penetration boundary conditions at 
the two ends require that the velocity vanish at x=0 and 
/. 


Following Ilinksii et al. 2 , expressing the variables in 
finite Fourier series in time, the equations (1) and (2) 
can be reduced to a system of Ordinary Differential 
Equations (ODEs) for the Fourier coefficients of the 
velocity potential cp, defined as u = V(p : 
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where 


v = (r 2 - r 2 )^,cp = • 

a = ^jd k e' h and 


a u = a ki (y k 'Vv *). fk = fk (v* .&.*)• The 
detailed expressions for a kl and f k are the same as those 
in Ilinksii et al. 2 except that the central blockage is 
taken into account. 
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where O , C 0 is the small signal 

Po 

propagation speed, D' m is given by 

K = ~(r - - im{y - 1 )(p m 
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and 

(v ! ],= itw,. 

l=-N+k 

The no-slip boundary conditions at the two ends are 
translated to the Fourier space: v k = 0 at x=0,l. 
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After the Fourier coefficients { (p k } and { v k } are 
obtained by solving the boundary-value problem Eq.(4), 
the velocity potential (p and the modified velocity v are 
computed from the inverse FFT. The density p is given 
by the momentum equation 

Po 


1 - 


Y - 1 


<v + l 

dt 2(r 2 -r 0 2 ) 2 


-v +ax -- 


S dv 
r 2 - r 0 2 dx 




Y - 1 


and the pressure p can be obtained from the state 
equation p! p 0 =(p/ p 0 ) r . 


II. Numerical Methods 


In this section, we describe the procedures for finding 
the optimal shape parameters so that the desired 
standing waveform can be obtained. 

To simplify our discussion, we introduce the 
following dimensionless variables 



T = a>t, R = j, R o =0j., 



where co is the frequency of the periodic force acted on 


7TC 

the resonator, coq = is the fundamental frequency 

of a cylindrical resonator of length / . 

In this work, we assume the acceleration of the 
resonator is harmonic, 

A(T) — Acos(iT) . The acoustic wave field in the 
resonator is determined by the acceleration amplitude 
A , the ratio of specific heats y, the attenuation 

coefficient G, defined as G = + , and the 

c oPo 


resonator oscillating frequency Q, defined as Q = 

co Q 

Now, we briefly describe the method for finding the 
optimal resonator shape. 

In the optimization process, we fix the values of the 
acceleration amplitude A = 5 x 10" 4 , the gas specific 
heat ratio y=1.2 and the viscosity-related parameter 
G=0.01. Suppose the resonator shape R(X) is 
determined by a number of shape parameters, S 0 , Sj, 

S n For example, a cone is given by 
R(X ) = S 0 + SjX . The compression ratio of a 

resonator defined at the ratio of maximum pressure 
to the minimum pressure at the narrow end of the 
resonator 


R c (S 0 ,S l S„)=-^ atX=0, (6) 

P min 

is a function of the shape parameters, the dimensionless 
frequency Q and the history of Q (due to the existence 
of hysteresis effects). The method for obtaining for 
fixed resonator shape is explained later in the section. 

We first outline the optimization schemes. A quasi- 
Newton method, BFGS (Broyden-Fletcher-Goldfarb- 
Shanno 4 ), is used for maximizing the multi-variable 

nonlinear function R c (S 0 ,S h S n ) . Since the 

evaluation of the objective function Rc itself involves 
solving a nonlinear system of ODEs Eq. (4) for many 
times, the gradient information of R* required for the 
BFGS method is not available analytically and is 
derived by partial derivatives using a numerical 
differentiation method via finite differences. This 
entails perturbing each design variables, S h in turn and 
calculating the rate of change in the objective 

function R c (S 0 ,S l S n ) . For two shape parameters, 

the optimization takes 4 to 48 hours of CPU time 
(depending on the type of the resonator shape and 
initial guess for the shape parameters) on a 1.3GHz 
Athlon T-Bird PC with the Lahey-Fujitsu FORTRAN 
compiler. 

For a given shape of the resonator and the central 
blockage, the boundary value problem Eq. (4) is solved 
numerically by a Multiple Shooting Method. Because 
of the hysteresis effects, the solution is not unique near 
the resonant frequency and the Multiple Shooting 
Method will not converge unless a good initial guess of 
the solution is provided. To circumvent the difficulty, a 
continuation method is implemented: the system of 
ODEs is solved starting from a frequency D. that is far 
away from the resonant frequency and the solution is 
used as an initial guess for solving the ODEs for 
increased or decreased Q; the steps are repeated until 
the all branches of the solution for all values of Q near 
the resonance is completed. The maximum ratio R c 
among different values of Q is chosen as the 
compression ratio for the resonator. 

III. Results 

A. Characteristics of the standing waves 

Before we present our results on finding the optimal 
resonator’s shape, we illustrate some of the important 
properties of the pressure wave in a non-cylindrical 
resonator, whose shape is described by that of a 
homcone: 
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[S 0 cosh(5,X),0< X <S 2 
\a + j3X,S 2 <X <1 * ' * 

where 

5 0 =0.028333, 5, = 5.7264, 5 2 = 0.25,a = 5 0 cosh(5,5 2 ), 
and p - 5 0 5, sinh^^). The shape parameters for 

the homcone are obtained from those used in the 
experiments by Lawrenson et al . 1 

In Fig. 2, we show the pressure waveforms at the 
ends of cylindrical resonator and the homcone 
resonator at their corresponding resonance frequencies 
for the effective viscosity G=0.01, the gas specific heat 

ratio y=1.2, the acceleration A = 10” 3 . Throughout 
this paper, we will keep the effective viscosity G and 
the gas specific constant at these values, which are 
close to the values that used in the experiments 1 and our 
numerical simulations have suggested that the pressure 
waveform is not a strong function of these parameters. 
At the same size of the acceleration, the waveform at 
the narrow end of the homcone shows large variation in 
which the ratio of the maximum and minimum 

pressure, R c , exceeds the value of 30.5. On the other 

hand, the waveform for the cylinder shows the 
formation of shocks at the resonance and the 
compression ratio R c is below 1.32. The difference in 

the waveforms shows the strong dependence of the 
acoustic field on the resonator shapes. The graph also 
shows that the variation in pressure at the wide end of 
the homcone is much milder than that of the narrow 
end, oscillating within 23% above or below the value of 

the reference pressure p 0 . 

In Fig.3, we plot the ratio of the amplitudes of the 
second and the first harmonics of the pressure, 

p 2 / P\> against the reduced amplitude of the first 

harmonic, p { / p 0 , measured at one end of the cylinder 

and at the narrow end of the homcone. The figure 
shows that the second harmonic reaches its maximum at 
small amplitude of the first harmonic for the consonant 
(cylinder) resonator and the second harmonic increases 
slowly as the first harmonic increases for the dissonant 
(homcone) resonator. In other words, the energy is 
pumped into the first harmonic more for the homcone 
than that for the cylinder, creating the large-amplitude 
pressure wave in the homcone resonator. In contrast, 
the energy is efficiently absorbed by all harmonics of 
the pressure wave in a cylindrical resonator and the 
amplitudes of the harmonics reach their maximums at 
small amplitude of the pressure. 

In Fig.4, the amplitude of the first harmonic is plotted 
against the frequency of the oscillating resonator for 


two different levels of acceleration: A = 2x10 4 and 

A = 1CT 3 . For the smaller acceleration, the pressure is 
uniquely determined at each frequency; for the larger 
drive amplitude, the pressure takes one of the multiple 
values near the resonance, depending on the direction 
of the change in frequency. The existence of hysteresis 
and hardening resonance in homcone resonators 
requires that the largest pressure compression ratio be 
obtained through an upward frequency sweep. Our 
numerical simulations show that the frequency 
increment size near the resonance must be small for the 
convergence of the Multiple Shooting method. 

B. Optimal conical shapes 

In following subsections, we present results on 
optimizing resonator shapes to achieve maximum 
compression ratio R c at one end of the resonators. For 
simplicity and reducing computing cost, we assume the 
central blockage is absent, i.e., fl 0 =0, and keep ratio of 
specific heats y=1.2, the attenuation coefficient G=0.01 
and the acceleration A = 5xlCT 4 constant. 

A conical resonator can be written as 
R(X) = S 0 +S x X , for 0 < X < 1 . (8) 

Using the optimization procedure described in Sec.II 
and starting with the conical shape given in Lawrenson 
et al. \ we found that the compression ratio R c of a 
conical resonator reaches a maximum when the cone 
has the shape constants: 5^0.032945 and 5i=0.26800. 
The achieved compression ratio is R c =5.048 for the 
fixed values of the parameters given earlier. Recall that, 
due to hysteresis, in order to find R Ci one has to trace 
out the entire branch by incrementing the frequency Q. 
The compression ratio reaches the value when the 
frequency is increased to 1.3135. The conical resonator 
studied in Lawrenson et al . 1 and Uinksii et al. 2 has 
almost identical shape 5^=0.032941 and 5 1 =0.26800. 

The adopted optimization scheme BFGS is designed 
for finding a local extreme of a multivariable function. 
Our numerical simulations indicate that the 
compression ratio, as a function of the shape 
parameters, have many local extrema. Starting with a 
different initial shape, we obtained another locally 
optimized conical shape with 5^=0.0104 13 and 
5j=0. 19572. When the frequency is raised to the value 
of 1.3862 by small increments, the compression ratio 
reaches the value of 5.343, which is about 6% higher 
than that for the first optimal conical resonator. In our 
optimization, we restrict the slope of the cone to be less 
than 0.268 so that the one-dimensional model used in 
this work is still valid. And, we require that the radius 
at the small end be larger than 0.01 so that real 
resonators can be built using these optimized 
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dimensions. The two optimal conical resonators are 
shown in Fig. 5. The second cone (Fig. 5(b)) has the 
smaller narrow end and smaller slope than the first one 
(Fig. 5(a)). 

C. Optimal horn-cone shapes 

Horn-cone geometry is described by 

Jtf.cosh^X) for 0<X<X c 
\a+j3{X-X c ) for X C <X<1 

where 

a = R h cosh(D h X c ) and /? = R h D h sinh (D h X c ) 

. R h , X c and D h are positive constants. We consider that 
the connection point of horn and cone in the homcone 
is fixed at X c =0.25. We optimize the remaining two 

shape parameters S 0 = R h and S x = D h . We found 

that the dimensions of the horn-cone given in 
Lawrenson et al.\ R h =0.028333, D h =5.7264, are 
nearly locally optimized if we start the optimization 
with these values. For these parameters, the maximum 
compression ratio 13.10 is achieved when the frequency 
is increased to 1 .467 1 . In searching the optimal 
homcone design, the parameters, such as acceleration 
and specific heat ratio, are fixed at the values given in 
the beginning of Sec. Ill B. Starting with a different 
initial geometry of the homcone, we can obtain a higher 
value of the compression ratio R c =17.12 when R h 
=0.024488, D h =5.1434. For the compression ratio, the 
frequency was increased to 1.41 10 gradually with small 
increments. Comparing with the compression ratio for 
the homcone in Lawrenson et al. x (shown in Fig.6(a)), 
the second optimal homcone (shown in Fig. 6(b)) 
improved the ratio by more than 30% at the same value 
of acceleration. Shown in Fig.6(c), the pressure 
waveform at the narrow end for new homcone design is 
a little more complicated than that for the original 
horncone given Lawrenson et al . 1 

D. Optimal V^-cosine shape 

We define a !/ 2 -cosine shape as 

R(X) = S 0 +S ] ( 1-cosOdO, for 0 < X < 1 . 

The parameters of the resulting optimal cosine shape 
are 5^=0.016307, S^O.093041, shown in Fig. 7(b), 
compared with So=0.025, Si=0.095, the ^-cosine 
dimensions reported in Chun and Kim 3 . For the optimal 
shape, we obtain the compression ratio Rc =10.67 at the 
frequency D=1.5683, which is about 11% better than 
that of Chun and Kim 3 at the same level of acceleration. 
Under the specified conditions, we found that hom- 
cone shape is better than the cosine shape in generating 
higher compression ratio at the narrow end, which is 
opposite to the findings of Chun and Kim 3 . Shown in 


Fig.7(c), the pressure waveform of the optimal */ 2 - 
cosine design is smoother than that of Chun and Kim 3 . 

E. The effect of the central blockage 

In this section, we discuss the influence of the 
central blockage on the pressure waveforms in a 
resonator of fixed shape. In this work, we limit 
ourselves to cylindrical central blockage, though our 
equations derived in Sec. I apply to arbitrary 
axisymmetric shape of central blockage. 

Figure 8 shows the pressure waves at the narrow end 
of a conical resonator described by Eq. (8) with 
So= 0.032945 and S^O.26800 for different central 
blockage radius sizes: R 0 = 0, 0.01 and 0.02. These 
pressure profiles were obtained at the fixed frequency 

fi= 1.3132 and for A = 5X10" 4 , G = 0.01, y = 1.2. 
The central blockage in this case was a simple cylinder. 
It appears that the presence of a central blockage 
reduces the amplitude of the pressure wave. However, 
Fig. 9 shows the resonant frequency of the conical 
resonator with a central blockage is shifted. Keeping 
other parameters constant, the pressure waveforms for 
the conical resonator with the cylindrical blockage of 
radius Rtf= 0.02, at different oscillating frequencies Q = 
1.2812, 1.3068 and 1.3291, are shown in Fig. 9. 
Comparing to the pressure wave shown in Fig. 8 for a 
resonator without a central blockage, the amplitude of 
the pressure is nearly the same when the central 
blockage of radius 0.02 is included. The results indicate 
that a cylindrical central blockage in a conical resonator 
would affect the resonant frequency but would have 
negligible effect on the amplitude and overall shape of 
the pressure wave at the resonant frequency of the 
resonator/blockage system. 

IV. Conclusion 

We have performed local optimization schemes for 
finding the resonator shapes that maximize the pressure 
compression ratio at one end of the resonators. The 
dimensions for several types of optimal resonators are 
reported: cone, homcone and ^-cosine shapes. For each 
type, we found there are many different designs that 
achieve local extrema. Trying with different initial 
guesses for the optimal design, we show that one can 
get as much as 30% improvement on the compression 
ratio with a fixed level of acceleration. For the shapes 
we considered, it appears that the homcone shape 
generates the highest compression ratio. Strategies for 
searching globally optimal shapes are under 
investigation. 

The effect of including a cylindrical central blockage 
in a resonator is studied in this work as well. For a 
conical resonator, it seems that the resonance frequency 
is increased when a central blockage is present and the 
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compression ratio at resonance is not significantly 
affected. 
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Figure 1 . Sketch of an axisymmetric resonator and a central blockage. 



Figure 2. The reduced pressure p / p 0 at the ends of a homcone and a cylinder at their corresponding 
resonance frequencies for the same acceleration A = 10~ 3 . 
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Figure 3. The ratio between the amplitudes of the second and the first harmonics of the pressure wave, 
p 2 l p x , is plotted as a function of the dimensionless amplitude of the first harmonic reduced pressure 
p x ! p 0 . These amplitudes are calculated for the narrow end of the homcone and at one end of the 

cylinder. 



1.35 1-4 1.45 1.5 


Q 

Figure 4. The effect of the frequency and the amplitude of acceleration on the dimensionless amplitude of 
the first harmonic reduced pressure p x / p 0 . The pressure is calculated for the narrow end of the homcone 

for the acceleration A = 2 X 10~ 4 (crosses) and A = 10~ 3 (circles and squares). 
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Figure 5. Optimized conical shapes: the compression ratios corresponding to the shapes in (a) and (b) 
reach the values of 5.048 and 5.343 respectively. For shape parameters, see the text. 


(a) (b) 



Figure 6. (a) The homcone shape in Lawrenson et al x (b) Optimized homcone shape. For shape 
parameters, see the text, (c) The pressure waveform at the narrow end for the homcone in Lawrenson et 
ai l (the solid line) and that for the second optimized homcone (the dashed line) are shown. 
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Figure 7. (a) The ^-cosine shape as in Chun and Kim 3 , (b) The optimized Vi-cosine shape, (c) The pressure 
waveform at the narrow end for the !/ 2 -cosine resonator in Chun and Kim 3 (the solid line) and that for the 
optimized ‘/ 2 -cosine shape (the dashed line) are shown. 



T 


Figure 8. Effect of the center blockage. The pressure waveforms at the narrow end of a cone resonator are 
plotted for the different radii of the central blockage: R 0 = 0 (solid line), ^=0.01 (dashed line) and 
i?o= 0.02 (dotted line). The parameters corresponding to the plot are given by Q = 1 .3132, A = 5 x 10^ , G = 

0.01, y= 1.2. 
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Figure 9. Pressure waveform at different frequencies: f! = 1.2812 (the dotted line), Q = 1.3068 (the dashed 
line) and f2 = 1.3291 (the solid line), fora conical resonator with the cylindrical central blockage of radius 

flo= 0.02. The rest of the parameters corresponding to the plot are given by A = 5 X 1 0" 4 , G = 0.0 1 , y = 
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C Software (Computer Codes) 

The computer codes for the project are included in the floppy disk attached 
to this final report. 
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